
	
	/*******************************************************************************
*   			Impacts of a Large-Scale Parenting Program: 				   *
*					Experimental Evidence from Chile						   *
********************************************************************************


	REQUIRES:	"${dt_final}panel_RR_indexes.dta"
	CREATES:	${out_tables}/T14A_Iv_parental_individual.xls
				
	WRITEN BY:  Italo Lopez  [italolop@usc.edu]

********************************************************************************
	Prepare data
*******************************************************************************/

	* Load data
	use "${dt_final}panel_RR_indexes.dta", clear


	* Establish macros for regressions

/*******************************************************************************
*   Set macros for regressions
*******************************************************************************/
	
	* Outcomes
	*-----------------------------------------------------------
	local childoutcomes  "zDCCSPunt zptevir_irtscore2 zt_score zcbcl_t_e zcbcl_t_i zibatt_sspi zibatt_ssai zibatt_sssr zas_raw_score" 
	local behaviors "home_score se_score disc_score" 
	local beliefs "zpstyle1 zpstyle2 zpstyle3 zPSCS zPSSS_family zPSSS_friends zPSSS_others zPSI_Distress zCESD zpacotis"	
	local child "zptevir_irtscore2 zcog_index zsoc_index" 

	* Controls
	*-----------------------------------------------------------

	macro  def Xvar1 "i.age_year i.gender"  	
	macro  def Xvar2 "i.age_year i.gender inc_qtaut_old i.pc_edu_mdsfin hh_mem hh_tipo"  	
	macro  def Xvar2old "i.gender inc_qtaut_old i.pc_edu_mdsfin_old hh_mem hh_tipo"  

	
	* Parental characteristics: WAIS and BFI and imputed wais 
	*-----------------------------------------------------------

	macro def parent_cha4 "BFI_ext BFI_ope BFI_con BFI_neu BFI_agr wais2 flag_wais"
	


/*******************************************************************************
	Table A14: IV estimates of parental individual measures
*******************************************************************************/
	
	local parental "zhome_score_gsem znur_score_gsem zdisc_score_gsem1b zdisc_score_gsem2b zPSCS zpacotis zPSSS_family zPSSS_friends zPSSS_others zIAP_authoritarian zIAP_authoritative zIAP_permissive zPSI_Distress zCESD" 
	#delimit;
	foreach var of local parental {;
				
					
				* PANEL 3: ITT with Age and Gender and SES;
					xi: xtivreg `var' $Xvar1 (treat2 treat3 = group2 group3) if x==1,  fe i(CENTRO_SALUD);
					estimates store `var'_iv1;
					
				* PANEL 4: ITT plus parental characteristics -; 	  
					xi: xtivreg `var' $Xvar2 $parent_cha4 (treat2 treat3 = group2 group3) if x==1, fe i(CENTRO_SALUD);
					estimates store `var'_iv2;
					
	};


	#delimit ;
	

	esttab zhome_score_gsem_iv1 znur_score_gsem_iv1 					  
		 zdisc_score_gsem1b_iv1 zdisc_score_gsem2b_iv1 zPSCS_iv1 		  
		 zpacotis_iv1 zPSSS_family_iv1 zPSSS_friends_iv1 				  
		 zPSSS_others_iv1 zIAP_authoritarian_iv1 zIAP_authoritative_iv1   
		 zIAP_permissive_iv1 zPSI_Distress_iv1 zCESD_iv1  				  
		 using "${out_tables}/appendix/T14A_iv_parental_individual.csv", replace f		   
		 keep(treat2 treat3) 											 
		 cells("b(fmt(3)star)" "se(fmt(3)par)") 							 
		 star(* 0.10 ** 0.05 *** 0.01) 								     
		 stats(r2_o N, fmt(2 0) labels("R-squared" "Observations")) legend	 
		 mtitles("Home Index" "PBC Nurturing Index"						  
		 "PBC Negative Discipline Index" "PBC Positive Discipline Index" 	 
		 "Perceived Self-efficacy" "PACOTIS scale" 	 
		 "Perceived Family Support" "Perceived Friends Support" 			  
		 "Perceived Support from Others" "Authoritative Style"			  
		 "Authoritarian Style" "Permissive Style" "Parental Mental Health" 
		 "Parental Stress""Depression" )  								 
		 title("PANEL 1: Iv with Age and Gender and SES") 	varlabels(treat2 "NEP-B" treat3 "NEP-I") ;
		

		esttab zhome_score_gsem_iv2 znur_score_gsem_iv2 					 
		 zdisc_score_gsem1b_iv2 zdisc_score_gsem2b_iv2 zPSCS_iv2  		  
		 zpacotis_iv2 zPSSS_family_iv2 zPSSS_friends_iv2 				  
		 zPSSS_others_iv2 zIAP_authoritarian_iv2 zIAP_authoritative_iv2 
		 zIAP_permissive_iv2 zPSI_Distress_iv2 zCESD_iv2				  
		 using "${out_tables}/appendix/T14A_iv_parental_individual.csv", append f		   
		 keep(treat2 treat3) 											  
		 cells("b(fmt(3)star)" "se(fmt(3)par)") 							  
		 star(* 0.10 ** 0.05 *** 0.01) 									  
		 stats(r2_o N, fmt(2 0) labels("R-squared" "Observations")) legend	  
		 mtitles("Home Index" "PBC Nurturing Index"						  
		 "PBC Negative Discipline Index" "PBC Positive Discipline Index" 	 
		 "Perceived Self-efficacy" "PACOTIS scale" 	 
		 "Perceived Family Support" "Perceived Friends Support" 			  
		 "Perceived Support from Others" "Authoritative Style"			  
		 "Authoritarian Style" "Permissive Style" "Parental Mental Health" 
		 "Parental Stress""Depression" )  								 
		 title(" PANEL 2: iv plus parental characteristics") 	
		 varlabels(treat2 "NEP-B" treat3 "NEP-I");
